Structural relaxations in electronically excited poly(para-phenylene) 
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Structural relaxations in electronically excited poly(para-phenylene) are studied using many- 
body perturbation theory and density-functional-theory methods. A sophisticated description of 
the electron-hole interaction is required to describe the energies of the excitonic states, but we show 
that the structural rela:xations associated with exciton formation can be obtained quite accurately 
within a constrained density-functional-theory approach. We find that the structural relaxations in 
the low-energy excitonic states extend over about 8 monomers, leading to an energy reduction of 
0.22 eV and a Stokes shift of 0.40 eV. 

PACS numbers: 71.15Mb, 71.15.Qe, 71.20.Rv, 71.35.Aa 
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The sustained growth of interest in conjugated poly- 
mers arises from their potential as the active material in 
low-cost field-efFect transistors, photovoltaic devices and 
light emitting diodes (LEDs). 1, 2] Investigating electron- 
ically excited states is an important part of this effort, 
which has proved a great challenge for theoretical tech- 
niques. An accurate description of the excited electronic 
states of conjugated polymers requires the inclusion of 
the strong electron-hole interaction and the structural 
relaxations which occur in the excited state. 

In its simplest form, a polymer LED consists of a layer 
of conjugated polymer sandwiched between two elec- 
trodes. Electrons and holes are injected into the polymer 
where they are attracted to one another and form bound 
excitons. Light is then emitted by exciton recombination. 
In an excited state the polymer may lower its energy by 
structural relaxations so that the peak optical emission 
frequency is lower than the peak absorption frequency. 
Although this "Stokes shift" can, at least in principle, 
be measured, the actual relaxations have not been de- 
termined experimentally. These relaxations localize ex- 
citons and are relevant for the technologically important 
processes of exciton migration and recombination. 

We have investigated structural relaxations in the low- 
energy excited states of poly(para-phenylene) (PPP). A 
PPP LED has been demonstrated which emits blue light 
in a band around 2.7 eV.0,0| PPP is, however, insolu- 
ble and therefore difficult to process and, instead, soluble 
derivatives of PPP with various side groups are preferred 
for manufacturing LEDs. The low-energy optical proper- 
ties of PPP and its derivatives are similar, although the 
structural relaxations in the excited states depend on the 
nature of the side groups. We have chosen to study PPP 
because it has a simpler structure than its derivatives. 

As we show below, systems containing more than 100 
atoms are required for studying the excited-state relax- 
ations of PPP, which pose a significant problem to theory. 





FIG. 1: The structure of a PPP chain with the C atoms 
represented by filled circles and the H atoms by open circles, 
(a) Side view showing the A and B bonds, (b) End view 
showing the torsional angle 9. 



Such systems can easily be treated by density-functional 
theory (DFT) 0, but DFT alone often has difficulties 
in describing excited states. Excited states can be de- 
scribed by many-body perturbation theory (MBPT) j^, 
but this is numerically demanding and cannot currently 
be applied to such large systems. However, we find that 
the forces in the excited states of PPP obtained within 
DFT are very similar to those within MBPT, which in- 
dicates that the relaxed structures and Stokes shift ob- 
tained within DFT are reliable. 

A PPP chain (Fig. ^) consists of phenyl rings joined 
by single C-C bonds, with adjacent rings being rotated 
with respect to one another by a torsional angle 9. Our 
DFT structures for the ground states of the single chain 
and crystal are in excellent agreement with a previous 
DFT calculation and in accord with experiment. Our 
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DFT band gap for the single chain is 2.44 eV and that for 
the crystal is 1.70 eV, which is in good agreement with 
a previous DFT calculation. 

We began by using a simple model of how an opti- 
cal excitation affects the bonding in PPP: we performed 
DFT calculations of the changes in the electronic charge 
distribution due to a single excitation from the highest 
occupied molecular orbital (HOMO) to the lowest unoc- 
cupied molecular orbital (LUMO). In both the isolated 
chain and crystalline forms we found an increase in the 
electronic charge in the A bonds and to a lesser extent in 
the B bonds (see Fig.P). These changes indicate the de- 
velopment of double-bonding character which is expected 
to favor both a more planar arrangement of rings (reduc- 
tion in 9) and compression of the A and B bonds. 

In a perfect chain or crystal an exciton would, in prin- 
ciple, be delocalized because localization within some 
region of space costs kinetic energy. However, a larger 
energy may be gained by structural relaxations so that 
an exciton may localize itself. The low-energy exci- 
tonic states of PPP may therefore be self-localized by 
structural relaxations, i.e., local reductions in 9 and the 
lengths of the A and B bonds. 

The structural relaxations in the lowest excited states 
of a single PPP chain can be investigated using a con- 
strained DFT method. To this end we model the excited 
state by promoting an electron from the HOMO to the 
LUMO at wave vector A; = and performing full struc- 
tural relaxations in unit cells of various sizes. A unit cell 
containing 28 monomers was found to be sufficient to 
converge the excited state structure. The important pa- 
rameters of the minimum energy structure are shown in 
Fig. [3 The relaxations are largely confined to a region of 
approximately 8 monomers over which the torsional angle 
and the lengths of the A bonds are significantly reduced, 
and the lengths of the B bonds are somewhat reduced. 
The gain in energy from the relaxation is 0.22 eV, and 
the Stokes shift, calculated as the difference between the 
band gaps of the ground and excited state structures, is 
0.40 eV.|i 

We have analyzed the localization of the exciton in 
terms of the compression of the A and B bonds and 
the reduction in the torsional angle, 9. We performed 
two constrained DFT calculations for the excited state 
in which (1) only relaxations along the axis of the chain 
and (2) only relaxations in the plane perpendicular to 
the chain were allowed. The relaxation energies of 0.14 
eV and 0.17 eV, respectively, add up to significantly more 
than the full relaxation energy of 0.22 eV, indicating that 
the modes are strongly coupled. The relaxations along 
the chain led to a Stokes shift of only 0.06 eV and to 
very weak localization of the exciton, whereas the relax- 
ations perpendicular to the chain gave a Stokes shift of 
0.16 eV and much stronger localization. The torsional 
relaxation is therefore very important in determining the 
Stokes shift, although the coupling to the bond compres- 



sion is strong and both types of relaxation are important 
in determining the full relaxation energy and Stokes shift. 
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FIG. 2: The torsional angle, 6, and the lengths of the A and 
B bonds (see Fig.0 as a function of monomer number for the 
excited state of a single PPP chain. 

Stokes shifts in PPP derivatives have been measured 
in dilute solutions and solid films. A common approach 
is to use "laddered" PPP derivatives in which some or all 
of the phenyl rings are joined to their neighbors by chem- 
ical "bridges". These bridges hinder rotation about the 
A bonds which, according to our picture, is expected to 
reduce the Stokes shift. Hertel et al. |^ report absorption 
and emission measurements on a series of laddered and 
non-laddered polymers, which clearly show that ladder- 
ing reduces the Stokes shift. As shown in Fig. 3b of Her- 
tel et al. Q , in dilute solution the non-laddered polymer 
dodecyloxy-poly(para-phenylene) (DDO-PPP) exhibits a 
Stokes shift of 0.6±0.2 eV. Considering the presence of 
the solvent and different side chains in the experimental 
system the agreement with our single chain PPP value 
of 0.40 eV is reasonable. The measured Stokes shifts for 
the laddered polymers are smaller, in agreement with our 
picture that changes in 9 couple strongly to the band gap. 
Our picture also suggests that structural relaxations for 
excitons in laddered PPP derivatives should only lead to 
a weak localization effect. 

We now investigate the reliability of the constrained- 
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DFT results for the relaxation and Stokes shift of the ex- 
citon. The formation of a completely delocalized exciton 
leads to a vanishingly small change in the charge den- 
sity and therefore vanishingly small forces on the atoms. 
If, however, a lattice distortion occurs which reduces the 
band gap in some region of space then the exciton will 
be attracted to this region. If the lowering of the energy 
of the exciton is greater than the energy required to form 
the distortion then the exciton will stabilize the distor- 
tion, leading to a self-localized exciton. The reduction in 
the band gap due to a lattice distortion is an effect which 
is described approximately within our DFT calculations. 
The excitonic energies themselves, however, are strongly 
modified by many-body effects arising from the electron- 
hole interaction. We now investigate the dependence of 
these many-body effects on the lattice distortion. If this 
dependence is sufficiently weak the constrained DFT ap- 
proach will be reliable. 

A rigorous approach to electronic many-body effects 
for excited electronic states is given by many-body per- 
turbation theory (MBPT).|3| We have used the GW ap- 
proximation to describe the addition or removal of an 
electron and the Bethe-Salpeter equation (BSE) for the 
excitation of an electron including the electron-hole inter- 
action. These techniques have rece ntly been used to de- 
scribe the excitonic states of solids [lOj . clusters ^3 1 ^^'^ 
polymers. MBPT is numerically very demanding 

and cannot currently be applied directly to a polymer 
with a self-localized exciton. 

To investigate the effects of electronic correlation and 
validate the constrained-DFT approach we proceeded as 
follows. We performed DFT and MBPT calculations for 
a series of structures involving torsional angles and com- 
pressions of the A and B bonds which are similar to those 
found in the self-localized exciton of Fig. |21 We consid- 
ered the relaxation of the two central monomers of Fig. [3 
as indicating the maximum relaxations which occur in 
the exciton. We then periodically repeated this struc- 
ture, obtaining a reference system which can be inves- 
tigated within both constrained DFT and MBPT. We 
defined a structural parameter x which takes the value 
x=0 for the ground state structure (0=33.7°, A=1.456 A, 
and B= 1.380 A), and x=l for the geometry correspond- 
ing to the central monomers of the self-localized exciton 
{0=9.1°, A=1.427 A, and B=1.366 A). Intermediate val- 
ues of X correspond to linearly interpolating 9, A and B 
between the extremal values. Fig.Olshows the excitation 
energies of this periodic structure as a function of x. The 
four curves denote the DFT energy gap, the quasiparticle 
(QP) energy gap, and the transition energies of the low- 
est spin-singlet and spin-triplet excitons obtained from 
solving the BSE.[l3 The QP gaps are 2.3-3.0 eV larger 
than the DFT gaps due to the significant QP corrections 
typical of semiconducting systems. The exciton energies, 
on the other hand, are 1.5-2.9 eV smaller than the QP 
gaps due to the attractive electron-hole interaction. The 



singlet (triplet) excitation energy ranges from 3.41 eV 
(2.62 eV) at a;=0 to 2.18 eV (1.61 eV) at x=l. 

The change in the excitation energy when going from 
x=0 to x=l constitutes the excited-state contribution to 
the force on the structure projected onto the relaxations 
studied. The most important feature of Fig.|31is that this 
change (-1.22 eV or -1.01 eV, respectively, for the singlet 
and triplet state) is very similar to the change in the DFT 
gap energy (-0.95 eV). Changes in the DFT gap energy, 
in turn, correspond to the excited-state forces given by 
constrained DFT. Our results therefore imply that the 
forces are indeed given reliably by the constrained DFT 
theory, even though the exciton binding energies are large 
and the absolute values of the DFT band gaps are in poor 
agreement with the excitonic energies. [15j 

The singlet excitation energy for the ground state 
structure of the single chain (3.41 eV) is in good agree- 
ment with experimental absorption spectra which show 
maximum absorption around 3.5 eV.0,0| We find that 
the singlet-triplet splitting depends rather weakly on the 
relaxations, changing from 0.79 eV at x=Q to 0.57 eV 
at x=l. We therefore expect the splitting for the re- 
laxed excited state structure to be in the region of 0.5 
eV, which is consistent with the splitting of 0.7-0.8 eV 
obtained from emission spectra of PPP derivatives, jsj 
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FIG. 3: Excitation energies calculated within the DFT and 
GW approaches and the GVK-BSE singlet and triplet energies 
as a function of the structural parameter x for a single PPP 
chain. 

We have analyzed the reasons for the success of the 
constrained DFT approach for calculating excited state 
geometries. The details of this analysis will be published 
elsewhere, but we conclude that the various approxima- 
tions should work well for the chain or solid when the 
excitons are delocalized over many atoms and when the 
relaxations are relatively small and do not include re- 
constructive changes such as bond breaking. The con- 
strained DFT approach would be inaccurate if the exci- 
tonic wave function contained contributions from more 
than one electron or hole band, but this mixing could 
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be calculated from BSE calculations on high symmetry 
structures and then used in DFT studies of large struc- 
tures. Our study suggests that the constrained DFT ap- 
proach is likely to work for other conjugated polymers. 

Our main calculations are for single chains, in which 
the electron- hole interaction is larger than in solids. 
Therefore our calculations are an even stronger test of the 
idea of neglecting the electron-hole interaction when cal- 
culating the excited state relaxations than would be en- 
countered in the more technologically relevant solid state. 
In the solid, many chains are packed together at van der 
Waals distances, which gives rise to three main additional 
features: (i) electronic overlap between chains, (ii) van 
der Waals interactions, and, in particular, (iii) the effect 
of inter-chain dielectric screening on the electron-hole in- 
teraction. The first of these is well described within DFT 
while MBPT includes the second and third. 

To study the effects of intcr-chain dielectric screening 
we considered a 3D crystalline array of PPP chains in 
which the electronic overlap between chains is small and 
consequently the DFT results are only weakly perturbed. 
The calculated QP correction to the DFT gap and the 
electron-hole interaction are, however, both reduced by 
more than 1 eV as a consequence of the more effective di- 
electric screening in a 3D solid. The QP state, which de- 
scribes an additional electron or hole on one chain, shows 
significant interaction with the polarizable neighboring 
chains, and the QP gap closes. In addition, the QP cor- 
rections also depend on x, indicating that the intra-chain 
screening changes when the geometry relaxes. The ex- 
citon energies are, however, only weakly affected by the 
changes in the inter- or intra-chain screening. In the solid 
the singlet (triplet) energy ranges from 3.66 eV (3.14 eV) 
at a;=0 to 2.55 eV (2.23 eV) at x=l, which is only a lit- 
tle higher than the single-chain results. This insensitivity 
derives from the charge-neutral character of the exciton, 
which is much less influenced by electrostatic screening 
than charged single-particle excitations. Similar effects 
have also been observed for other polymers by van der 
Horst et al. 0] The changes in the excitation energies as 
X varies from to 1 (-1.11 eV for the singlet, -0.91 eV for 
the triplet) are very similar to the change of the DFT gap 
(-0.95 eV), confirming the reliability of the constrained- 
DFT forces. The weak interactions between chains have 
to be handled with care, however, and the self-trapping 
of the exciton in the solid state will be the subject of a 
further study. 

In summary, we have shown that while the electron- 
hole interaction is very important in determining the ex- 
citonic energies in PPP it is less important for the excited 
state geometries. We have demonstrated that a simple 
constrained DFT approach is sufficient for calculating re- 
laxations in the low-energy excited states of PPP. In this 
case we find that the structural relaxations extend over 
about 8 monomers, leading to an energy reduction of 
0.22 eV and a Stokes shift of 0.40 eV. The DFT method 



is tractable for studying the influence of structural disor- 
der on the optical properties of polymers. 
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